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In multidimensional gravity with an arbitrary number of internal Ricci-flat factor spaces, interacting with electric 
and magnetic p-branes, spherically symmetric configurations are considered. It is shown that all single-brane black 
hole solutions are stable under spherically symmetric perturbations, whereas similar solutions possessing naked 
singularities turn out to be catastrophically unstable. The black hole stability conclusion is extended to some classes 
of configurations with intersecting branes. These results do not depend on the particular composition of the D- 
dimensional space-time, on the number of dilatonic scalar fields ip a and on the values of their coupling constants 
P_h ' Asa. Some examples from 11-dimensional supergravity are considered. 

m 

(N 

1. Introduction 

O , 

In this paper we continue our studies of multidimensional gravitational models based on D -dimensional Einstein 
equations with fields of antisymmetric forms of arbitrary rank (see [1-3] and references therein) as some low-energy 
limit of a future unified model (M-, F- or other type). Our main interest here will be in the stability properties 
of multidimensional black-hole (BH) and non-BH solutions with nonzero fields of forms, associated with charged 
p-branes. There exist a large number of such solutions in arbitrary dimensions — see e.g. [2, 4-9] and references 
therein. They are important in connection with studies of processes at early stages of the Universe, counts of 
micro-states in BH thermodynamics and now especially due to new developments in M-theory ]10[ related to the 
O !* AdS/CFT correspondence p| . For recent reviews of this rapidly developing field see, e.g., []l2|, [l3| . 

BH stability studies have a long history, of which we will only mention (more or less arbitrarily) some mile- 
stones, concerning spherically symmetric backgrounds. Regge and Wheeler |lj] considered the stability of the 
Schwarzschild space-time and developed the formalism of spherical harmonics for metric perturbations. Vishvesh- 
wara [jl5| finally proved the linear stability of Schwarzschild BHs; Moncrief |ll| did the same for Reissner-Nordstrom 
ones. BHs with a conformally coupled scalar field were shown to be unstable under spherically symmetric per- 
turbations fll7||, as well as minimally coupled scalar field configurations in general relativity possessing naked 
singularities ]18[. The monopole degree of freedom is present there due to the scalar field; it was argued that 
monopole perturbations were most likely to be unstable due to the absence of centrifugal terms in the effective 
potentials; catastrophic instabilities were indeed found and it was unnecessary to study other multipoles. On the 
other hand, coloured BHs, containing non-Abelian gauge fields, were shown to be, in general, unstable due to their 
sphaleronic degrees of freedom — see |l9| and references therein. A recent overview of 4-dimensional perturbation 
studies may be found in Ref. |^0| . 

For BHs in multidimensional theories of gravity the situation is more complex since, on the one hand, there 
emerge new effective scalar fields (extra-dimension scale factors, sometimes called moduli fields) in the external 
space-time, and, on the other, instabilities may be caused by waves in extra dimensions. Instabilities of the latter 
kind were indeed found by Gregory and Laflamme [ pTj , p2| for a limited class of neutral and charged black strings 
and branes, having a constant internal space scale factor. Furthermore, it was argued that compactification on a 
sufficiently small length scale should prevent the onset of instability, and, moreover, that extremal black branes 
are stable ^3|. It was concluded that only very light BHs, whose horizon radii have the same order of magnitude 
as their extra dimensions, manifest this form of instability. 

It is therefore of interest to inquire whether or not there are other forms of instability, maybe "more dangerous" , 
on more general backgrounds, containing nontrivial internal space structures and/or several dilatonic scalars and 
brane charges. As was previously the case with backgrounds containing effective scalar fields, it is natural to 
consider first the simplest, monopole perturbations. 
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Earlier we analyzed the stability of static, spherically symmetric solutions to the Einstein-Maxwell-scalar 
equations with a dilatonic type coupling between scalar and electromagnetic fields in D -dimensional gravity p4j , 
p5| . It was proved there that only BH configurations were stable under linear spherically symmetric perturbations, 
while non-BH solutions turned out to be catastrophically unstable. A similar result was obtained for dilatonic 
BHs with the inclusion of the Gauss-Bonnet curvature term due to one- loop quantum corrections |2(| . We will 
now show that in the simplest case of a single charged black brane the solution is stable under linear spherically 
symmetric perturbations, whereas single-brane solutions with naked singularities are unstable. So the results of 
P4 p5[ are generalized. 

We also present a tentative consideration of multi-brane BHs and conclude that in cases when the perturbation 
equations decouple, the stability conclusion is also valid. Two classes of such systems are indicated, both charac- 
terized by certain relations among brane charges, such that, in terms of Sec. 3, the constituent vectors Y s form a 
single block of a block-orthogonal system (BOS) — single-block BHs for short. Namely, the stability is proved for 
arbitrary two-brane single-block BHs and multi-brane single-block BHs with mutually orthogonal vectors Y s (see 
the details in Sec. 6.2). For many single-block configurations which do not belong to these classes, the st abili ty 
can be proved as well, but their properties require individual studies; see an example in the Appendix, Eqs. ( A. 4 )- 
(A.7). There are, however, numerous multi-brane BHs for which decoupling is impossible and one may expect 
that some of them show a new type of instability connected with mode interaction; a study of these systems is in 
progress. 

The paper is organized as follows. Sec. 2 describes the general features of the field model to be considered. 
Sec. 3 presents some known static solutions, including BHs, on the basis of the target space V connected with 
dimensional reduction. In Sec. 4 a truncated target space V, more appropriate for treating the perturbations, is 
introduced, and wave equations for perturbations are derived. In Sec. 5 the stability properties of single-brane 
configurations are deduced, while in Sec. 6 the stability of some multi-brane BHs under spherically symmetric 
perturbations is established. The Appendix gives some examples from 11-dimensional supergravity. 

The word "stable" throughout the paper means "stable under linear spherically symmetric perturbations" . 



2. The model 

Our starting point is, as in Refs. [1-8], the model action for D -dimensional gravity with several scalar dilatonic 
fields tp a and antisymmetric n s -forms F s : 

S = J d D Zy/W\{m- 5a b9 MN d M ^d N V b - Y,^* 2 ^^}, (!) 

in a pseudo-Riemannian manifold M = R„ x Mo X ... x M„ , with factor space dimensions di, i = 0, . . . ,n; 1Z 
is the scalar curvature. We will assume M to be spherically symmetric, so that the metric is 

n 

ds 2 D = g MN dz M dz N = e 2a °du 2 + £ e 2 ^ ds 2 

n 

= c 2a °du 2 + c 2 ^°dn 2 - e^dt 2 + ]T e 2 ^ds 2 . (2) 

Here u is a radial coordinate ranging in R„ C R; ds^ = dfl 2 is the metric on a unit do -dimensional sphere 
Mo = S ; t G Mi = M t is time; the metrics g % = ds 2 of the "extra" factor spaces (i > 2) are assumed 
to be Ricci-flat and can have arbitrary signatures £j = signg 1 ; |<?| = |detgA/./v| and similarly for subspaces; 
F 2 = F s . Mi...M n Fs 1 " s ; A sa are coupling constants; r] s — ±1 (to be specified later); s S S, a € A, where S 
and A are some finite sets. The "scale factors" e' 3 and the scalars ip a are assumed to depend on u and t only. 

The _F-forms should be also compatible with spherical symmetry. A given F-form may have several essentially 
(non-permutatively) different components; such a situation is sometimes called "composite p-branes"^. For conve- 
nience, we will nevertheless treat essentially different components of the same F-form as individual ("elementary") 
F-forms. A reformulation to the composite ansatz, if needed, is straightforward. 



3 There is an exception: two components, having only one noncoinciding index, cannot coexist since in this case there emerge 
nonzero off-block-diagonal components of the energy-momentiim tensor (EMT) , while the Einstein tensor in the l.h.s. of the 
Einstein equations is block-diagonal. See more details in Rcf. [fll . 
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Each n s -form F = dA = 9[Af 1 -^-M 2 ...M n ]dz Ml . . . dz Mns is then associated with a certain subset / = . . . , i^} 
(i\ < . . . < ik) of the set of numbers labelling the factor spaces: {i} = Iq = {0, ...,n}. The forms F s are 
naturally classified as electric (F c j) and magnetic (F m i) ones. By definition, the potential Aj of an electric form 
-F c / carries the coordinate indices of the subspaces M i; i d I and is u -dependent (since only a radial component 
of the field may be nonzero). A magnetic form F m j is built as a form dual to a possible electric one, and its 

nonzero components carry coordinate indices of the subspaces M.;, i £ I = Iq \I , One can write: 

n e j = rankFe/ = d{I) + 1, n m / = rank F m j = D — rank .F e j = d{I) (3) 

where d(I) = J^iei are the dimensions of the subspaces M/ = X . . . X M{ fc . The index s will be used to 
jointly describe the two types of forms, so that Q] 

5 = {s} = {e7 s }u{mJ s }. (4) 

We will make some more assumptions to assure that all F-forms behave like genuine electric or magnetic fields 
in the physical subspace M p hy S = K« x Rt x Mo, namely: 

(i) 1 e I s , Vs (the subspaces M Is contain the time axis R t ); (5) 

(ii) ^ I s , Vs (the branes only "live" in extra dimensions); (6) 

(iii) —Tt(F a ) > 0, Vs (the energy density is positive). (7) 

By (i), the so-called quasiscalar forms (2| |j (forms with f I s , behaving as effective scalar or pseudoscalar 
fields in M p h ys ) are excluded. The reason for adopting (i) is that our interest here is mostly in BHs which do not 
admit nonzero quasiscalar forms (the no-hair theorem for brane systems ||). 

Assumption (iii) holds if all extra dimensions are spacelike (sj = 1 , i > 2 ) and in (Q) all r) 8 = 1 . In more 
general models, with arbitrary Ej , (iii) holds if 



Vd s - -e(Is), Vmi s = e(Is), e{I) = JJ £i . (8) 



def 

iei 



We will consider static configurations and their small (linear) time-dependent perturbations. It turns out, 
however, that under the above assumptions the Maxwell-like field equations for F s may be integrated in a general 
form for their arbitrary dependence on u and t. Indeed, for an electric m-form F s (s = el , m — d(I s ) + 1) the 
field equations due to (|l|) 

are easily integrated to give 

F utM 3 ...M m ^ Qse - Q °~* -2X saV ° £ M 3 ...M d(I) , A^T ^ l~F^ = £(I)Q 2 s e- 2 <^- 2X '^ a . (10) 

* to! 

where £ " and e... are Levi-Civita symbols, \gi\ = Yiipj Id^l' an< ^ Q s = cons t are charges. In a similar way, for a 
magnetic m-form F s (s = nxT, m = d(I s )), the field equations and the Bianchi identities dF s — lead to 

F s>Ml ...M da) = Qse Ml ...M d(1) J\^\ => ^F* = e(l)Qle-^ +2X ^\ (11) 



s,Mi...M d(7 j = Vs£jlfi...M (j ( T) Y \yj\ => —^s 

We use the notations 

n 

cn = Y J d 3 P°(u,t), a(I) =^d i /3 i (u,t). (12) 

j=i i£l 

Evidently, the expressions ([!(]) and ( |Tl| ) differ only in the signs before A sa and the signature-dependent pref- 
actors e. Due to (||), their energy- momentum tensors (EMTs) coincide up to the replacement A sa — > — A sa , and 
their further treatment is quite identical. In what follows we therefore mostly speak of electric forms, but the 
results are easily reformulated for any sets of electric and magnetic forms. We also assume that all Q s ^ 0. 
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3. Static systems 
3.1. The target space V 

Under the above assumptions, the system is well described using the so-called a model representation p]), to 
be briefly outlined here as applied to static, spherically symmetric systems. This formulation can be derived by 
reducing the action (|f|) to the (d + 1) -dimensional space R„ x M . 

As in [^7j and many later papers, we choose the harmonic u coordinate (V m Vm« = 0), such that 

n 

oP{u)=a Q {u) = Y, d ^- ( 13 ) 
i=0 

Due to (||), the combination (\) + ( 2 ) of the Einstein equations has a Liouville form and is integrated giving 

!kr x sinh ku, h > 0, 
u, h = 0, (14) 

fc _1 sin ku, h < 0. 

where k is an integration constant. With ( |l4| ) the D -dimensional line element may be written in the form 

(d d = f d -i) 

-2o-i/rf f j 2 "I " 



[ds(k,u)] 2 / d I [ds(fc,u)] 



i=2 



The u coordinate is defined for < u < u max where u = corresponds to spatial infinity while u max may be 
finite or infinite depending on the form of a particular solution. 

The remaining set of unknowns /3 % (u), ip a {u) (i = 1, . . . ,n, a 6 A) can be treated as a real-valued vector 
function x A (u) (so that {^4} = {1, . . . , n} U A) in an (n + \A\) -dimensional vector space V (target space). The 
field equations for x A can be derived from the Toda-like Lagrangian 

n 2 

l = c AB xt< + E ^ {U) = E d ^uf + + + E e 2Vs(u) (16) 

s i=l s 

(the subscript u means d/du), with the "energy" constraint 

E = G AB x A u x B u -Y,Ql e 2 ^ = -r^V 2 sign k. (17) 
^ — ' do — i 

s 

The nondegenerate symmetric matrix 
{GAB)= ( d ^rd + d i 5 ij \ 



<U, 

defines a positive-definite metric in V; the functions y s {u) are scalar products: 

y s = a(I s ) - \ sa ip a = Y StA x A , (Y S>A ) = fan,, -A so ), (19) 

where 5u — 1 if i £ / and <5;/ = otherwise. The contravariant components and scalar products of the vectors 
Y s are found using the matrix G AB inverse to G A b ■ 

ft/di - 1/(D - 2) \ /v M _ f. d{I s ) 



s y.y ,, = n I,,) - + (21) 

The equations of motion in terms of Y s read 

x A = J2Q 2 Ys A e 2VsH - (22) 
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3.2. Exact solutions: orthogonal systems (OS) 

The integrability of the Toda-like system (jl^) depends on the set of vectors Y s . In many cases general or special 
solutions to Eqs. ( p^ ) are known. Here we will mention the simplest case of integrability: a general solution is 
available if all Y s are mutually orthogonal in V 0, that is, 

Y S Y S , = 5 SS ,Y 2 , Y 2 - d(I) [1 - d(I)/(D - 2)] + A 2 > (23) 

where A 2 = ^ A 2 a . Then the functions y s (u) obey the decoupled Liouville equations y s ,uu — Q 2 S Y 2 e 2ys , whence 

e -2y s (u) =QlY 2 S 2 (h s , U + U S ) (24) 

where h s and u s are integration constants and the function s(., .) has been defined in For the sought 

functions x A (u) and the "conserved energy" E we then obtain: 

s s 

E = j : !^h 1+52= do k2signk! (26) 

where the vectors of integration constants c and c are orthogonal to all Y s : c A Y s a — c A Y s a = 0, or 

c l d t S lIs - c a \ sa = 0, c 1 d l 5u s - c a \ sa = 0. (27) 

3.3. Exact solutions: block-orthogonal systems (BOS) 

The above OS solutions are general for input parameters (D, di, Y s ) satisfying Eq. (p3|): there is an independent 
charge attached to each (elementary) F-form. One can, however, obtain special solutions for more general sets of 
input parameters, under less restrictive conditions than (^3|). Namely, assuming that some of the functions y s (u) 
( |l9[) coincide, one obtains the so-called BOS solutions jlj, where the number of independent charges coincides with 
the number of different functions y s (u) . 

Indeed, suppose Q that the set S splits into several non-intersecting non-empty subsets, 

S = |JS W , \S u \ = m(w), (28) 

UJ 

such that the vectors Y n(u) (m( w ) G Sw) form mutually orthogonal subspaces V w C V: 

Y Ku) Y v(un =0, lo + J. (29) 

Then the corresponding result from Q can be formulated as follows: 

Proposition 1. Let, for each fixed ui, all Y v e N u be linearly independent, and let there be a vector Y u = 
S/js5„ a M^M with a M > such that 

YvYu = Y 2 = Yl , V/i e 5 W . (30) 



Then one has the following solution to the equations of motion (|22|), ([I 

xA = E^-( w ) + cAu +^' ( 31 ) 

e" 2 ^ = q u Y 2 s 2 (h u ,u + u u ), q u d = Y] Q 2 (32) 



yr^hlsignhu 2 d 2 . 
E = 2^ y% ^do~l k Slgnk (33) 

where h u , u u , c A and c A are integration constants; c A and c A are constrained by the orthogonality relations 
(^) (so that the vectors c and c are orthogonal to each individual vector Y s 6 V). 

Eqs. ( |30| ) form a set of linear algebraic equations with respect to the "charge factors" a v = Q 2 /q u , satisfying 
the condition Ylues a M ~ 1- A solution to ( |j0|) for given Y M can contain some a M < 0; according to Q], this 
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would mean that such a p-brane is "quasiscalar" , violating the assumption (|5j). Solutions with such branes are 
possible but are rejected here since they do not lead to black holes. Furthermore, if a solution to (|3(]) gives 0^ = 
for some this means that the block cannot contain such a p-brane, and then the consideration may be 

repeated without it.Q 

The function y u (u) is equal to y^^ju) — Y^^A xA , which is, due to (|30|), the same for all [i 6 S u . The BOS 
solution generalizes the OS one, (p4|), (p5|): the latter is restored when each block contains a single J 7, -form. 

Both kinds of solutions are asymptotically flat, and it is natural to normalize the functions y s (u) and y w (u) 
by the condition y s (0) = or y^(0) — 0, so that the constants u s and u w are directly related to the charges. 

Other solutions to the equations of motion are known, connected with Toda chains and Lie algebras ^, @, |9) ■ 



3.4. Black-hole solutions 

Black holes (BHs) are distinguished among other spherically symmetric solutions by the existence of horizons 
instead of singularities in the physical do + 2 -dimension space M p h ys ; the extra dimensions and scalar fields are 
also required to be well-behaved on the horizon to provide regularity of the D -dimensional metric. Thus BHs are 
described by the above solutions under certain constraints upon the integration constants. Namely, for the BOS 
solution (|30|)-([33|), requiring that all the scale factors e' 3 * (except e' 9 = -\/l3tt| which should tend to zero) and 
scalars tp a tend to finite limits as in w max , we get Q : 

h u = k>0, Vw; c A = kJ2 Y J 2 ^ A ~ kSt (34) 



where A = 1 corresponds to i = 1 (time) . The constraint (|26|) then holds automatically. The value u = w max = oo 
corresponds to the horizon. The same condition for the OS solution (p^)-(p7|) is obtained by replacing whs. 
Under the asymptotic conditions ip a — > , (3 l — > as u — > , after the transformation 

e -2fc» = 1 _J* d d ^d -l (35) 

dr d 

the metric ( |l5| ) for BHs and the corresponding scalar fields may be written as 

„2 



ds 



i=2 



. def 2 ^ a M d(J M ) os 2 <f(J g ) 4 def 2 r os 2 

A " = ^2 E = y^^-^ ; A " = ~y2 E a ^ = -y2^; (36) 



' a = -E^ ln ^ E «a» °= s -Eyf 1 ^- ( 37 ) 

where == means "equal for OS, with lo hh s", and iT^ are harmonic functions in R + x S d ° : 



fl*(r) = 1 + ?W(dr d ), Pw d = v^ 2 + 4^ - fc- (38) 



The subfamily (g4J), (|36[)-(|3§|) exhausts all BOS BH solutions; the OS ones are obtained in the special case of 
each block consisting of a single element s . 

The above relations describe the so-called non-extremal BHs. Extremal ones, corresponding to minimum black 
hole mass for given charges (the so-called BPS limit), are obtained in the limit k — > 0. The same solutions follow 
directly from (|3^)-(|32]) under the conditions h u = k = c A = 0. For k = 0, the solution is defined in the whole 
range r > , while r = in many cases corresponds to a naked singularity rather than an event horizon, so that 
we no more deal with a black hole. However, in many other important cases r = is an event horizon of extremal 
Reissner-Nordstrom type, with an AdS near-horizon geometry; some examples are mentioned in the Appendix. 

Other families of solutions, mentioned at the end of the previouis section, also contain BH subfamilies. The 
most general BH solutions are considered in Ref. Q. 



4 Geometrically, the vector Y ^ solving Eqs. ( ^o| ) is the altitude of the pyramid formed by the vectors Y M , £ Suj with a common 
origin. The condition > means that this altitude is located inside the pyramid, while a M = means that the altitude belongs to 
one of its faces. 
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4. Perturbation equations 
4.1. Truncated target space V 

Consider now nonstatic spherically symmetric configurations corresponding to the action ([!]) with the metric (Q) 
and all field variables depending on u and t . As before, we are dealing with true electric and magnetic forms F s , 
so that their I s 9 1 , or 

7 S = 1UJ S) J s C{2,...,n}. (39) 

As in Refs. [^4| p5| , it is helpful to pass to the Einstein frame in the physical (d + 2) -dimensional space-time 
Mphys = R u X Mo X Mi . The action ([!]) is then rewritten in terms of the metric g^ u , the do + 2 -dimensional 
part of gMN , and is transformed to the Einstein frame in M p h ys with the metric 

V = e 2 <^S„- (40) 
The electric n s -forms are re-parametrized as follows: 

F u tM 3 ...M ns = F ut , — F 2 = ii^F^e- 2 ^) = -Q 2 s e~ 2X ^\ (41) 

s n s \ I s s 

where the indices M3, . . . M„ 3 belong to J s ; here and henceforth the indices /i, v are raised and lowered using the 
metric <?„„; in the last equality the solution ( p!o| ) and the positive energy assumption (J?]) are taken into account. 
The action (|lf) is written in terms of ~g^ v and F ^ as follows (up to a constant prefactor, connected with the 

volume of extra dimensions, and a subtracted total divergence): 



S = 



r I n 

\g\i [ TZ[g]--(da 2 ) 2 -J2di(dl3 i f 



i=2 

5 ab (d<p a , d<p b ) - i y F^F^ e W*>-2»W+2W 

Z Z — ' S S 



f d d ° +2 z^W\{n[g]-H KL {dx K 1 dx L )- l -Y, F ^ Ve ~ 2Z 

J { Z ses s s 



1 phys 



(42) 



where {df,dg) = g^ v d^fd v g , (<9/) 2 = (df,df); the non-degenerate symmetric matrix 



(H KL ) (** 'V" A ' I ) l«) 



b 

defines a positive-definite metric in the vector space V (truncated target space) parametrized by the variables 
(x K ) — (/3 2 , . . . , (3 n ; ip a ) ; the constant vectors Z s S V are characterized by the components^] 

(Z S . K ) = (dtSij. - J, -A sa ), (Z S K ) = (H KL Z S>L ) = - -X sa ) (44) 

where the matrix (H KL ) is inverse to (Hxl)'- 

{HKL) = (sv di -i/(D-2) oy (45) 



r KM_ (8*/d i -\/{D-2) 

^ S ab y 

The truncated target space V may be considered as the hyperplane x 1 = — 02/do in V, with the metric Hkl 
induced by Gab (^D)- The components H KL turn out to be the same as G AB for i ^ 1; the components Y S A 
and Z S K coincide in the same manner. It is easy to find that for vectors whose I s satisfy ( |39] ) (which is always 
the case in the present paper), 

Y S Y S ,^Z S Z S , + '^-^-, (46) 
d 

whence it follows that, first, when different Y s are mutually orthogonal in V, the corresponding Z s are never 
mutually orthogonal in V; second, for any Y s whose I s 3 1 one has 

Y S Y S ^Y 2 >^-1. (47) 
do 



5 We will use the indices K,L for quantities specified in V to distinguish them from those in V where the indices A,B are used; 
vectors in V are marked with overbars, those in V by arrows. Scalar products are written as YZ = GabY a Z b (as before) and 
YZ = H KL Y K Z L . 
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4.2. Wave equations 

The action ( pt2| ) may be used to obtain the equations governing small spherically symmetric perturbations of static 
solutions. The metric ( p0| ) in M p h ys will be written in the form 

dsl = g^dz^dz" = e 2a du 2 + e 2p dfl 2 - e 2 ">dt 2 (48) 

where "E" stands for the Einstein frame and a(u,t), (3(u,t), j(u,t) are connected with the corresponding quan- 
tities from (j|) as follows: 

a = a°+a 2 /d , f3 = P° + a 2 /d , j = /3 1 +a 2 /d . (49) 

Since the field equations for the F -forms have been integrated — see ( |l0| ) and (^lj), the remaining unknowns are 
a,/3,7 and x (that is, (3 l , i>2, and tp a ). In what follows we will write 

a(u, t) = a (it) + 8a(u, t), 

where 8a is a small perturbation, and similarly for other unknowns. We accordingly preserve only terms linear in 
8a and similar quantities and in time derivatives. The field equations may be written in the form 

□ [9] x K = E Z - K Q> e- 2d0/3+2 ^, (50) 

K = t; -^t* t { f; (51) 

where □ [#] = g^V^Vv is the D'Alembert operator, while for the nonzero components of the EMT T£ corre- 
sponding to ( |42] ) one has (no summing in fi) 

^=c- 2 ^L«diag(0, 1, 0,...,0)-^Q 2 c- 2 «+ 2 ''+ 2 ^diagfl-l, -i, -1 ~), 

f ut = H KL x%xi (52) 

where x u = d u x and xt = dtx; the first and second places under the symbol "diag" belong to t and it, respectively. 
As in our previous papers on stability, we use the coordinate freedom in the perturbed space-time and put 

8(3 = (53) 



but preserve the harmonic it coordinate condition in the unperturbed (static) space-timeQ Then Eqs. ( |50| ) and 
© give 

L 8x K + x*{8 lu - 8a u ) - 2x*Ja = 2 £ Q 2 S c 2 ^ +2 ^Z s K Z. %L Sx L ; L = - e 2d ^d tt + d uu] (54) 

s 

d 8a t = x u 8x t ; (55) 
d f3u(8a u - Sju) = 2x uu 8x - 2f3 uu 8a, (56) 



where (|54| ) follows from the ( u j) component of ( pl| ) and fl56|) from one of the angular components of (|5l|); we have 
also used the equations valid for static systems, in particular, (p2|), where, according to the definitions of Y and 
Z , y s (u) = Y s x = 7 + Z s x. Integrating ( |55| ) in t and omitting the emerging arbitrary function of it (since we 
neglect static perturbations), we obtain 

do8a = x u 8x. (57) 

Substituting Sa from ( p7| ) and 8a u — 8^ u from ( |56| ) into ( p4[ ) , we finally arrive at the set of wave equations for 
the dynamical degrees of freedom in our system, represented by 8x K : 

L8x K = 2P K L 8x\ P K L = + y2Q 2 s e 2 y^ZfZ^ L . (58) 

The stability problem is now reduced to a boundary- value problem for 8x K (u,t). Namely, if there exists a 
nontrivial solution to Eqs. ( p8|) satisfying some physically reasonable conditions at the ends of the range of u, such 

"This coordinate is harmonic for both metrics (b|) and (0); in the latter the coordinate condition has the form a = do (3 + 7. 
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that \Sx K \ (at least some of them) grow unboundedly with t, then the static system is unstable. Otherwise it is 
stable in the linear approximation. 

The condition at infinity, u — 0, is evident: the perturbations must vanish, 

dx K as u — ► 0. (59) 



It is less evident at u = u max since some of the background static solutions are singular there. As in Refs. 18, 25]] 
and others, dealing with minimally coupled or dilatonic scalar fields, we will use the minimal requirement providing 
the validity of the perturbation scheme, namely 

\5x K /x K \<oo. (60) 

When the background is regular, this condition requires that the perturbation should be regular as well. 



5. Stability properties of single-brane solutions 



5.1. Decoupling cases 

Eqs. ( |58| ) in general do not decouple. Even in the simplest case when there is only one antisymmetric form F (that 
is, one p-brane), so that Eqs. (|24|)-(p7|) yield the general static solution to the field equations, Eqs. ([38]) contain 
various linear combinations of Sx K with it-dependent coefficients. 

There is, however, an important case when Eqs. (|5^) do decouple for any configuration of M with the metric 
(|[). namely, the single-brane solution (24)-(|2^) under the condition that the vector c = (c K ) is parallel to Z in 

B = const (61) 



BY K /Y 2 , 



(here and henceforth in this section we omit the index s since, by our assumption, it takes only one value). This 
condition is automatically valid for the case of utmost interest, BHs with one p-brane ("a black p-brane"), which, 
by j34|), corresponds to B = k = h > 0. 

Due to the collinearity condition ( |6l| ) and the constraint ( p6| ) , the constants are now connected by the relation 

N'(h 2 sign h - B 2 ) = k 2 sign k- B 2 (62) 



with N' = Y 2 (do — l)/do < 1 . It turns out that, besides BHs, the condition (61) is satisfied for some singular 
solutions whose behaviour is quite generic for the system under study: 

1. k > 0, h > 0, such that it max = oo and a singularity at the centre of symmetry is attractive at least in 
terms of the metric (fi8]) , e 7 — > ; 

2. h < 0, so that the solution behaviour is determined by the function s(h,u + ui) — h~ x sinh(u + ui) in 
( pi| ) where u\ = const € (0, ir/\h\). In this case the central singularity is repulsive, e 7 — > oo, of Reissner- 
Nordstrom type. 



Due to (|6l|), Eqs. ( |58| ) take the form 



L5x 



K 



2Z 



K 



j,, s dof y 
f(u) = 



1 t f 2 \ 1 - 

T ( I f) -fun (ZSX), 

with y(u) determined by (|2~4|); the area function j3 has the form 

13= __J_{ln[(do -l)s(k,u)] + —f(u)-Bu + const 
do - 1 I N' 



Bu 



Y 2 



(63) 



(64) 



where the value of the constant is inessential. 

Since V is an I -dimensional Euclidean space (I = n— 1 + \ A\ ), there are I — 1 linearly independent vectors 
Z±_ such that Z±Z = 0. Therefore the set of wave equations (|63|) decouples into one equation for Zbx and I — 1 
equations for different Z±Sx: 



L (Z5x) 
L (Zx&x) 



U{u)(Z5x), 
0. 



1 



U{u) = 2Z 2 — ^ -/„„ ; (65) 

(66) 

7 Curiously, they are parallel in V although the corresponding vectors c and Y are mutually orthogonal in the surrounding target 
space V . For a clear picture, imagine two vectors in 3-dimensional space whose projections onto a plane lie on the same ray. 
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Z ± Sx = ^'(u) e n \ 



The static nature of the background solution makes it possible to separate the variables: 
ZSx = ip{u)e n \ 
so that Eqs. ( |65| ) and (|6^) lead to 

t/w = [e 2d °"n 2 +U(u)]i/>, 



e^fi'V. 



(67) 

(68) 
(69) 



The existence of an admissible solution of any of these equations with a real value of Q or f2' would mean that 
the perturbation can grow exponentially with time, hence the instability. 

It is hard to solve Eqs. (68), (|69| ) in their full range but it is rather easy to assess the asymptotic behaviour of 
their solutions near u = and u max , and this will be sufficient for making stability conclusions. 

In particular, for u — > 0, which corresponds to spatial infinity, one has 



U{u) -> 



and 



e dof3 m C(s u- da ' d , 



CO 



jdo/d 



(70) 



The general asymptotic form of solutions to ( pq ) and (|6^) at small u for all cases under study may be written as 
follows: 



i/> (or V') = u do/{2d) 



c\exp(codVLu 1 ^ d ) + C2 exp(— codQu 



ci, C2 = const. 



(71) 



The boundary condition ( |59| ) then requires that in (71) c\ = 0, and it remains to look at the other end of the u 
range, u — > u max . 



5.2. Instability of naked singularities 

Consider Case 1 of the previous subsection, a "scalar type" singularity. As u — > oo , the relevant functions of the 
static solution behave as follows: 

y = hu + 0(1), (i « - ^^L_ l (fc _ _ B ) (72 ) 

where iVi = d^Y 2 /d > 1. Since, due to (|62|), in the present case 

\B\ > k <=^- \B\> h and \B\ < k <^ \B\ < h, 

one sees that — > exponentially. The same happens to {J (u) , therefore the asymptotic form of ( |68| ) or ( |69| ) 
is simply = or ip' uu — , so that 

(or ?/)') = c 3 u + c 4 (73) 

with constants C3 and C4 . On the other hand, the background functions x K also behave as const ■ u as u — > 00 , 
therefore the second boundary condition (^) is satisfied for any solution (|73|), including the one joining the solution 
( fFll) with ci = at small m . We conclude that there are growing modes of perturbations for any f2 , hence the 
singular solution is catastrophically unstable. 

In Case 2, h < 0, we have u max = 7r/\h\ — u\ < 00 and the relevant functions in the static solution approach 
u ma x in the following way: 

Z K 1 
y(u) w -]n(\h\Au), x K « -— \n{\h\Au), (3 w —— In Au -00 (74) 

where Au = u max — u. One can make sure that U(u) does not affect the asymptotic behaviour of solutions to 
Eq. ( |68| ) as u — > u max as compared with that of Eq. (|69|), and for both one can write: 

ip (or ^ ) = C5 + ceAu (75) 

while the condition ( |60| ) only requires Ifa^/ In Au| < 00. Thus the solution satisfies ( |60| ) for any choice of the 
constants C5,C6, and, as in Case 1, this leads to the instability of the background singular solution. 
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5.3. Stability of black holes 



In the BH case it is again hard to solve Eqs. (p8|), (|6S|), but for our purpose it is sufficient to note that, to realize 
an instability, a solution should begin with a zero value at it = and tend to a finite limit as u — > oo . This 
is evidently impossible for a solution to ([39]) since i/)' uu /ip' > 0. We conclude that at least the ip 1 modes of BH 
perturbations are stable. The same reasoning works for the ip mode provided U(u) > for all u > . Let us pass 
to the variable R — r d /d in the expression for U in (]65|), so that < u < oo corresponds to oo > R > 2k: 



U = 



2Z 2 pR(R - 2k) 
Y 2 (R + p) 2 
T i dd 1 d Q 



N 



2k- 
1 



1 



N'R' 



Y 2 d 



< 1, 



{R + P') 2 

, def 
P = 



pN' 
p(l-N% 



AkR + 2k(p + p') + pp' 
(R + p') 2 



(76) 



and the substitution ([55]) with 
, so that, in particular, p' > 0. 



where we have used the explicit form of the single-brane BH solution (J24|)- 
R = r d /d, replacing Y s — > Y , p s — > p, u s — > u\. Note that N' < 1 due to 

The expression ([76]) is manifestly positive for oo > R > 2k, therefore the ip mode also does not lead to an 
instability. Thus linear stability of all single-brane BH solutions under spherically symmetric perturbations has 
been established. 

Our consideration did not apply to extremal BHs since in this case the behaviour of the background functions 
x K (u) is generically singular as u — > oo (i2 = 1/u —> 0): 



(77) 



and so there is no reason to require \Sx K \ < oo. In some cases it is regular (see Sec. 3.2 and examples in the 
Appendix). One can see, however, that again, as u — > oo, Eqs. ( |68| ) and ([39]) for a single-brane extremal BH take 



the form ip uu — ; the linearly growing solution is discarded since it grows faster than x K in ( |77| ) , so we are left 
with a constant and have to require \if> < oo| for both regular and singular backgrounds. Then the same reasoning 
with ip U u/ip > makes us conclude that such allowed solutions with f2 > do not exist and extremal BHs are 
stable as well. Indeed, an explicit form of U (u) is 



U(u) 



2Z 2 



Y 2 (u + u 1 y 



1 + N 



,-u 2 + N"u 2 



(ill + N"u) 2 



(78) 



where N' < 1 was defined in (|7^) and N" = 1 — N' . The reasoning works since U > for all u > and U — > 

as u — y oo . 

We can now formulate the following result, to be used in the further consideration: 

Proposition 2. If a decoupled linear perturbation mode £ of a static, spherically symmetric BH solution obeys 
the equation 



with U{u) > (including the case U = ), this mode is stable. 



(79) 



6. Some black holes with multiple branes 
6.1. Two-brane black holes 

We have seen that one-brane singular background solutions are catastrophically unstable; we would not like to 
treat more complex singular solutions since there is no reason to believe that interaction of modes can prevent 
the instability. We instead consider some multi-brane BH solutions for which the perturbation equations decouple 
and show that they are stable. 



Suppose there is a BH background solution (p4[)-(38) with two branes, so that s takes two values, s = 1,2. 



The solution is characterized by two charges Q s , two vectors Y s e V and their counterparts Z s e V , which we 
assume to be non-collinear (if they are collinear, the consideration simplifies and the result is the same as for a 
single brane). 
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The matrix P K \ in the perturbation equations (|58| ) may be written in the form 
P K l =J2Q" e 2y ° {u) Z? Z S<L + £ ZfZ s ,, L fsAu), (80) 

3 SS 1 

with 



fas' (u) 



Y?Y?, 



(y s ,u + k)(y s ' <u + k) 



(81) 



Just as in the one-brane case, one easily separates the "transversal" degrees of freedom: for vectors Z± 6 V 
such that Zj_Z s — (they fill a (dimV— 2) -dimensional plane), the function £ = Z±Sx obeys the wave equation 
© with U = 0. 

However, Eqs. (|58|) with the matrix (]8(]) in general do not decouple. An exception is the special case when the 
two functions y s coincide, 

yi = y-2 = y(u) = ksmh{ku\) / sinh[fc(u + u\)] (82) 

although the vectors Y s are different; in the BOS terminology (Sec. 3.3) the two vectors Y „ form a block and 
in our case this single block exhausts the whole system. If we suppose for simplicity that the norms of Y s 
coincide, Y 2 = Y 2 2 — Y 2 , then the charges coincide as well, Q\ — Q\ = Q 2 , and one obtains for the two modes 
i± = (oZi ± Z 2 )5x: 

L£+ = 2Z 2 (l + cos6>) [Q 2 e 2y + 2F] C+; (83) 
= 2Z 2 Q 2 c 2y {\ -cos6»)£_ (84) 

where Z 2 = Y 2 — (d — l)/d (see (p6|), 8 is the angle between the vectors Z\ and Z 2 in V and F is the function 
@ with Y 2 = Y 2 and Vl = y 2 = y . 

A direct substitution of y and f3 u into ( |S3"1 ) shows, as before, that the coefficient by 8x + at the r.h.s. is 
nonnegative; for Sx~ in ( ]8"i| ) this is manifestly so. Hence the previous reasoning works and we conclude that such 
BHs (including extremal ones) with two branes are stable. 

The case Y 2 ^ Y 2 2 is covered in the next section. 

6.2. Single-block black holes 

A natural question arises, whether or not the stability conclusion of the previous section extends to an arbitrary 
multi-brane BH described by a single function y(u) , in other words, to any single-block BH. Note that any set of 
linearly independent vectors Y s may be treated as a BOS-block, hence a special static solution of this kind (and 
hence a BH solution) may always be obtained; the only restriction is a M > for the charge factors obeying the 
consistency conditions (|30"|). 

Consider such a system: let there be a BOS BH solution with m linearly independent vectors Y s G V, 
s € S = S u , and the charge factors a s satisfy (BOh. The following relations are valid: 



Y u =J2aaY s ; Y S Y^=Y 2 , Vs; = (85) 

s s 

It is easy to see that, due to (filf), similar relations hold for the corresponding vectors Z s 6 V: 

Z w =^a s Z s Vs; Z S Z U = Z 2 J . (86) 

s 

For certainty we suppose that Z s are linearly independent; if they are not, the consideration is slightly modified 
without changing the results. 

The wave equations (pq) take the form 



\L8x k = qe 2 y^Y, a sZ s K Z sL Sx L + F{u)Z^ K Z u L 5x L , F(u) - ^ 



(Vu + k) 2 



do ftu 



(87) 



where q = = J2 S Q 2 and y(u) is given by J8^). As before, the perturbations Zj_5x (where Z± belong to the 
plane Vj_ orthogonal to all Z s , whose dimension is dim V — m > ) are decoupled and obey the equation (|6^) , 
giving no unstable modes. 
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Multiplying ( p7|) by Z u , one obtains a decoupled equation for £ w = Z^Sx: 

= 2l7 w (u)^, I7 w (u) = (go 2 * + F)Z 2 

where = Z^,. Since, as is directly verified, U^(u) > 0, this mode is also stable. 

The remaining (m — 1) degrees of freedom may be described in terms of the vectors W s = Z s — Z^ and the 
functions £ s = W s 8x, such that 

W s Z UJ = 0; ^a s W s = 0; ^a s & = 0. (89) 

s s 

Using and (|39|), one obtains the following m equations, coupled due to (|S9|), for (m— 1) independent variables: 

m 

L £ s = 2q e 2 y £ , K ss , = ajW.W* ■ (90) 

s'=l 

Excluding one of the unknowns, say, £ TO , by virtue of (^), we arrive at a determined set of wave equations for 
Vs = 6 - £m , s = 1, . . . , m - 1: 



1 / i m - 1 \ 

Li h = 2qe 2y Y] F SS , V ' S , F ss , = a s ,W s , [W s + — V . (91) 

V a ™ s „ =1 / 



s' = l 



This is a good way of studying specific models. In the general case, however, the situation looks more transparent 



if we consider, instead, an auxiliary system with m independent unknowns, described by Eqs. (90) where W m is 
slightly shifted from its true value by some AW, so that all W s become linearly independent; the relation among 
l; s in ( ]89j ) is then cancelled as well. Our system is restored when AW — > . 

For the auxiliary system the matrix (W S W S >) is symmetric and positive-definite; if all a s are equal, the same is 
true for the matrix of coefficients in (pJ), (K ss /) , hence there is a similarity transformation bringing it to a diagonal 



form with its positive eigenvalues along the diagonal. Such a transformation applied to Eqs. (POJ) decouples them 



into m separate wave equations like (88), with some positive function replacing U(u). In the limit AW — > 0, 
the worst thing that can happen is that some of the eigenvalues tend to zero, giving for some combinations of £ s 
the equation L£ = which, as we know, does not lead to an instability. One can assert "by continuity" that 
this picture is generic, at least for a s close enough to being equal, and stability is again concluded according to 
Proposition 2. 

On the contrary, when the numbers a m > arc different, one cannot guarantee that the non-symmetric matrix 
K ss i is similar to a diagonal one p§| . A failure in its diagonalization can be connected with the occurrence of a 
pair (or pairs) of complex roots X s of the characteristic equation det \K SS > —XS SS > \ = 0. In this case there is at least 
one pair of coupled perturbations for which a special investigation is necessary. An inspection of the characteristic 
equation shows that the matrix K ss i cannot have negative eigenvalues, therefore a separate unstable mode cannot 
occur and the only possible instability can be connected with coupling between modes. 

In particular, in an arbitrary OS BH solution there is a subfamily where all y s {u) coincide (i.e., the constants 
m s are the same for all s), so that the branes form a BOS block, and it turns out that all a s are also equal, as 
well as the squared charges Q 2 S . The above reasoning shows that such solutions are stable. 

If rank(i^ ss /) < m — 1, that is, there are additional linear dependences among Z s , then some combinations of 
£ s decouple leading to equations of the form L £ = , and for the remaining modes the above discussion can be 
repeated with slight modifications. 

This is what can be said about the general case of single-block BOS BH solutions. If there is a block of only 
two branes (m = 2), one can make a common stability conclusion generalizing the one made in Sec. 6.1. Indeed, 
for £ w = + a 2 £,2 there is Eq. ({88|), whereas for = £i — £2 one obtains 

L£_ =2qe 2 y ai {l-a 1 )(Z 1 -Z 2 ) 2 ^. (92) 

In the special case Z\ = Z\ = Z 2 one recovers (p3|), (pi}). 
For m > 3 one has to study specific models individually. 
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7. Concluding remarks 

We have shown that all static single-brane BH solutions with the metric (||) are stable under linear spherically 
symmetric perturbations, whereas non-BH solutions possessing naked singularities of different types are unstable. 
Very probably other singular solutions, for which perturbation equations do not decouple, are unstable as well, 
since, as known from vibration theory, coupling between modes can hardly stabilize them. On the contrary, coupled 
modes can be unstable even when single ones are stable. It is therefore of interest to study the stability properties 
of more complex BH solutions; this work is in progress. 

We have also shown that the BH stability conclusion can be extended to some BHs with multiple intersecting 
p-branes, namely, for the BOS case, characterized by a single function y(u) (i.e., all y s coincide). It turns out that 
for such backgrounds the wave equations for perturbations also generically decouple and the absence of unstable 
modes can be proved. Though, such a general proof is available only in two cases: (i) two-brane BH solutions 
(to = 2) and (ii) equal-charge subfamilies of arbitrary OS solutions. Nontrivial brane systems with m > 2 should 
be studied individually to see whether or not the corresponding matrix (K ss i) in ( |90| ) can be diagonalized over 
the field of real numbers. If yes, the solution is stable, otherwise a further study of coupled modes is necessary. 

It should be stressed that a BOS-block solution exists for an arbitrary set of linearly independent vectors Y s . 
In particular, if any multi-brane static BH solution for a certain set of input parameters with independent vectors 
Y s is known, e.g., any OS or BOS solution (see Sec. 3), then the additional requirement that all the functions y s 
coincide selects from it a special BOS-block solution, for which a stability study can be performed as described 



above. The only restriction is the requirement a^{oj) > for the charge factors obeying Eqs. (30). 

Some technical points are worth mentioning. First, in gravitational stability studies it is sometimes rather 
hard to separate real physical perturbations from purely gauge degrees of freedom. We avoid this problem by 
obtaining the set of wave equations ( |58| ) where the number of equations is precisely the number of dynamical 
degrees of freedom, represented by scalars in the physical space-time M p h ys ■ Due to the latter circumstance, one 
more complication is avoided: when dealing with vector and tensor perturbations of BHs, one has to take into 
account the apparent singularity of the metric on the horizon; to properly formulate the boundary conditions, it 
is then necessary to pass to Kruskal-like coordinates; to be admissible, and the perturbations are required to be 
finite on the future horizon [|l5| ^l], p2[ . In our case the perturbations are scalars, so the finiteness requirement 
can be imposed in any coordinates. The choice of gauge only remains important for making the treatment more 
transparent. 

To conclude, we would like to emphasize that our consideration did not depend on the number and dimensions 
of the factor spaces in the original space-time M , on the number of scalar fields ip a and on the particular values 
of their coupling constants X sa ■ 



Appendix 

Consider, for illustration, some solutions of D = 11 supergravity, representing the low-energy limit of M-theory, 
as examples of systems to which our stability results apply. 

The action ([!]) for this theory does not contain scalar fields (<p a = X sa = 0) and the only F-form is of rank 4, 
whose various nontrivial components F s (elementary F -forms according to Sec. 2, to be called simply F -forms) 
are associated with electric 2-branes [for which d(I s ) = 3] and magnetic 5-branes [such that d(I s ) = 6] (see |lj] 
and references therein). The orthogonality conditions ( p3| ) are satisfied if the following intersection rules hold: 

303 = 1, 3H6 = 2, 6H6 = 4. (A.l) 

(the notations are evident); for all F-forms Y 2 = 2. 

We will designate the branes by figures labelling their world volume coordinates (covered by I s ) , beginning 
with "1" which corresponds to the time axis. Thus, e.g., (123) is an electic 2-brane whose world volume includes 
the time axis Mi = R t and two extra dimensions. The number of dimensions where branes can be located is 
D - 1 - d = 10 - d ■ 



1. Single-brane BH solutions are described by (pq), ( |3q ) where all sums and products in s consist of a single term. 
These solutions are well known; the metric (0) can be presented as 



ds 2 n = H d ^/ 9 



1 - 2k/{dr d ) 2 / dr 



H \l-2k/(dr d ) 



dt z + | — = -= r + rW ) + H-^ds^ + ds 2 oS 



(A.2) 
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where H = H(r) = 1 + p/(dr d ), p = yfk 2 + 2Q 2 — k, d = do — 1; ds 2 n and ds 2 s are the "on-brane" and 
"off-brane" extra-dimension line elements, respectively; the dimension do of the sphere Mo varies from 2 to 7 
for d(I) — 3 (an electric brane) and from 2 to 4 for d(I) = 6 (a magnetic brane). In particular, the cases of 
maximum do , when off-brane extra dimensions are absent, correspond in the extremal near- horizon limits to the 
famous structures AdS± x S 7 (electric) and AdS-j x S 4 (magnetic). All these solutions are stable under linear 
spherically symmetric perturbations. 

2. Some examples of orthogonal systems (OS), whose stability in the general case is yet to be studied, are: 

(i) (123), (145), (167) — 3 electric branes; d = 2 or 3. 

(ii) (123), (124567) — 1 electric and 1 magnetic branes; d — 2 or 3. 

(iii) (123), (145), (124678), (135678); d = 2. 

The metrics are easily found from ([56]) with D — 2 = 9 , Y 2 = 2 and the equalities marked ^ . 

The systems (i) with do — 3 and (iii) are remarkable in that their extremal limits have regular horizons and the 
near-horizon geometries are, respectively, AdS2 x S 3 xT e and AdS2 x S 2 x T 7 if the remaining extra dimensions 
are compactified on tori. 

When the orthogonal systems form BOS-blocks (i.e., in the special case of equal charges and a unique function 
y(u)), the solutions are stable according to Sec. 6.2. 

3. All two- brane BOS-block BHs are stable according to Sec. 6.1, for instance, 

(i) (123), (123456); Y X Y 2 = 1. (ii) (123), (145678); Y X Y 2 = -1. 

The norms are equal (Y 2 — 2), and the angle 9 is 60° in case (i) and 120° in case (ii). 

4. Many seemingly possible three-brane blocks turn out to be forbidden due to a zero value of a certain a s (see 
Sec. 3.3). Consider, e.g., 



(123), (145), (123678); (Y S Y S >) = 2 -1 . (A.3) 




One easily finds from (30) that (eti, a 2l as) = (0, 1/2, 1/2), so one of the charges should be zero, which means 
that such a system cannot exist. 

5. The following is an example of a single-block BH whose stability can be established by an individual study as 
described in Sec. 6.2: 




(123), (145), (123467); d = 2 or 3; (OV.! = | 2 " J ■ (A.J) 

From (^) it follows 

( ai , a 2 , o 3 ) = (2/7, 3/7, 2/7); Y u = ^a s f s ; Y 2 U = 6/7. (A.5) 



3 



s=l 



Suppose for certainty do = 3 . Then, in agreement with (Eq) 



4/3 -2/3 1/3 
-2/3 4/3 -2/3 
1/3 -2/3 4/3 



(Z...Z..m= [ -2/:-! 4/3 -2/3 | . Z'i = i-. (A.6) 



The next stage is to separate the perturbation Z^Sx, after which the remaining two degrees of freedom obey 
Eqs. @ of the form 

L Vs = 2ge 2 ^ £ F ss , Vs ,, F ss , = ( ^ . (A.7) 

The characteristic equation det \F SS > —X5 SS / 1 = has the form (2/7— A)(6/7 — A) = 0, and, according to Proposition 
2, the positivity of its roots proves the stability of the background configuration. 
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